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ABSTRACT 

We present a systematic search for disk features in 476 Virgo cluster early-type dwarf (dE) galaxies. 
This is the first such study of an almost-complete, statistically significant dE sample which includes 
all certain or possible cluster members with tob < 18 that are covered by the optical imaging data 
of the Sloan Digital Sky Survey Data Release 4. Disk features (spiral arms, edge-on disks, or bars) 
were identified by applying unsharp masks to a combined image from three bands (g, r, i), as well 
as by subtracting the axisymmetric light distribution of each galaxy from that image. 14 objects are 
unambiguous identifications of disks, 10 objects show 'probable disk' features, and 17 objects show 
'possible disk' features. The number fraction of these galaxies, for which we introduce the term dEdi, 
reaches more than 50% at the bright end of the dE population, and decreases to less than 5% for 
magnitudes tob > 16. Although part of this observed decline might be due to the lower signal-to- 
noise ratio at fainter magnitudes, we show that it cannot be caused solely by the limitations of our 
detection method. The luminosity function of our full dE sample can be explained by a superposition 
of dEdis and ordinary dEs, strongly suggesting that dEdis are a distinct type of galaxy. This is 
supported by the projected spatial distribution: dEdis show basically no clustering and roughly follow 
the spatial distribution of spirals and irregulars, whereas ordinary dEs are distributed similarly to 
the strongly clustered E/SO galaxies. While the flattening distribution of ordinary dEs is typical for 
spheroidal objects, the distribution of dEdis is significantly different and agrees with their being flat 
oblate objects. We therefore conclude that the dEdis are not spheroidal galaxies that just have an 
embedded disk component, but are instead a population of genuine disk galaxies. Several dEdis display 
well-defined spiral arms with grand design features that clearly differ from the flocculent, open arms 
typical for late-type spirals that have frequently been proposed as progenitors of early-type dwarfs. 
This raises the question of what process is able to create such spiral arms - with pitch angles like 
those of Sab/Sb galaxies - in bulgeless dwarf galaxies. 

Subject headings: galaxies: dwarf — galaxies: structure — galaxies: clusters: individual (Virgo) — 
galaxies: fundamental parameters — galaxies: spiral — techniques: image process- 
ing 



1. INTRODUCTION 

At first glance, early-type dwarf galaxies (dEs) are 
characterized by their smooth appearance, having no re- 
cent or ongoing star formation and apparently no gas 
or dust content. Since they are the most numerous 
type of galaxy in clusters, it is self-evident that most 
of the proposed formation scenarios for dEs reflect the 
vigorous gravitational forces acting within the very envi- 
ronment in which these gal axies typica lly reside. Ram- 
press ure stripping ( Gunn fc Gotd 119721) . g alaxv harass- 
ment IjMoore et al.ll996i) . and tidal stirring l|Maver et al.l 
|2001|1 are all based on the removal of gas and the morpho- 
logical transformation of a late-type spiral or irregular 
galaxy, thereby attempting to reproduce the seemingly 
plain appearance of dEs. On the other hand, differences 
in the chemical abundances of early-type and late-type 
galaxies m a,v argue against a s imple morphological trans- 
formation l)Grebel e t al''200^. In any case, such struc- 
tural transformations would be well-s uited to explai n the 
famous morphology-density relation IjDressleil 19801) : the 
higher the density is, the more efficiently are infalling spi- 
rals and irregulars transformed into dEs, thereby skew- 
ing the relative abundance of different types of galaxy 
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towards massive early-type objec ts as compared to abun- 
dances in the field. Moreover, IConselice et alJ l)200lD 
point out that the number of Virgo cluster dEs is more 
than a factor of 3 larger than what would be expected 
from just adding groups to the cluster. This strongly 
favours the idea that the majority of dEs were formed 
through a morphological transformation of galaxies that 
fell into the cluster. 

Especially in recent years, small or intermediate-sized 
samples of early-t ype dwarfs have been studied in a large 
variety of ways. iBoselli el al. (2005) find the relation 
of far-UV-near-UV colour and luminosity to behave op- 
posite for early-type dwarfs and giants. Van Zee et al. 
(2004a) derive intermediate ages and subsolar to solar 
metallicities for dEs via optical multiband photometry. 
Similar values were reported bv Geha et al. (2003) from 
a Lick index analysis of high-res olution spectra . These 
spectra and similar studies by Ivan Zee et alJ (|2004hD 
and Simien & Prugniel (2003) ^-Iso revealed a significant 
amoun t of rotation in some dEs. Finally, iBuvle et alJ 
()2005l) presented HI 21cm line observations as a first 
study of the interstellar medium of a dE outside the Lo- 
cal Group. 

However, no formation scenario could yet be clearly 
confirmed or rejected. This might be due to a very ba- 
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sic piece of the puzzle still lacking: the unambiguous 
characterization of early-type dwarf morphology. Fol- 
lowing common definition, early-type dwarfs comprise 
both dwarf ellipticals and dwarf SO (dSO) galaxies - we 
are not considering the fainter dwarf spheroidal galax- 
ies (ejj;. iGrebel et al. 2003) or the ultra-compact dwarfs 
(e.g.|Hilkcr ct al. 1999) here. The morphological appear- 
ance and overall profile of a dwarf elliptical are clearly 
defined. In contrast, dSO galaxies are loosely defined as 
objects whose overall appearance is similar to that of 
a dwarf elliptical, but where a more detailed examina- 
tion shows non-elliptical properties such as lens shap e or 
(central) asymmetries. iBinggeli fc CameronI (|199lD ar- 
gued that most of these characteristics were indicative of 
a disk nature, and the authors conjectured that "many, if 
not most, dSO systems must be disk galaxies" . However, 
their existence as a separate class of objects has been pu t 
in question by several authors fe.g. iRvden et al.lll999D . 
and dSOs have frequ ently been treated a s a subclass of 
dwarf ellipticals fe.g. 'Baraz za et al]l2003D . 

The unambiguous discovery of disk substructure (spi- 
ral a rms and/or bars) in some dwarf elli pticals and 
dSOsJ.Terien et al.ll2C)0ntl Baraz za et al.ll20n2l:lGeha et all 
l2003HGraham et al.ll200 3: De R iicke et al.ll2003(l event~ 
ally proved the presence of a disk in at least some early- 
type dwarfs. At the same time, however, this raised the 
question of whether these objects are genuine disk galax- 
ies, i.e. of flat oblate shape and without significant stellar 
spheroid, or whether they are spheroids hosting just a 
small disk component like the two low-luminosity ellipti- 
cals pres ented bylMorcUi et a l. (20 04). On the theoret- 
ical side, iMastropietro et all l)2005j) showed that a frac- 
tion of the progenitor galaxy's disk is able to survive the 
morphological transformation from galaxy harassment, 
providing a possible explanation for disks in early-type 
dwarfs. 

Since up to now, a systematic analysis of a large sam- 
ple of early-type dwarfs for the presence of disk features 
has been lacking, common practice has been to continue 
using the original classification of the Virgo cluster cat- 
alog (VCC, Binggeli ct al. 1985), therefore calling some 
objects 'dwarf elliptical', some 'dSO', and some 'dwarf el- 
liptical with embedded disk'. In order to avoid confusion 
we assign the common abbreviation 'dE' to early-type 
dwarfs in general, thereby meaning both dwarf ellipticals 
and dSOs. We shall then examine each object for poten- 
tial disk substructure, and introduce the term 'dEdi' for 
a dE with disk features. 

Clearly, the small sample of dEdis discovered so far 
can neither serve as basis for a revised classification nor 
is it sufficient to feed formation theories with quanti- 
tative input concerning the fraction and properties of 
such objects. A systematic search for disk features 
in dEs is thus required, and is made possible by the 
Sloan Digital S kv Survev (SP SS) Data Release 4 (DR4, 
lAdclman-McC axthv et aLl l2006'l which covers almost the 
whole Virgo cluster with multiband optical imaging. 
With these data at hand, our study can properly address 
the following questions: a) whether all objects listed as 
dSO in the VCC indeed show disk features, b) how large 
the fraction of galaxies with disk features is among dEs, 
c) how this fraction is distributed with respect to lumi- 
nosity, d) where in the cluster these objects are located, 
and e) whether they appear to be genuine disk galaxies, 



or just spheroids with a disk component. The catalog of 
dEdis and dEdi candidates resulting from this study will 
serve as important input for all future work on dEs, since 
the observables under study (e.g. dE colours) can then 
be correlat ed with the presence or absence of a disk. 

Recently, lAguerri et al.l l)200fit) have introduced a two- 
component definition of a dSO based on one-dimensional 
profile fits, with those (Coma cluster) objects being 
called dSOs where a single Sersic fit did not lead to a 
satisfying result and instead a combined Sersic plus expo- 
nential fit was necessary. Our goal in this paper, in con- 
trast, is to uncover disk features on the two-dimensional 
image without any presumption on one-dimensional pro- 
file shapes. To investigate whether or not the two defini- 
tions go hand in hand is beyond the scope of this paper, 
since it requires that accurate profile fits be done for all 
our SDSS galaxies. This will be the subject of a future 
paper in this series. 

Our data and sample selection is described in Sect. El 
followed by an outline of the techniques for image anal- 
ysis in Sect. |31 Identifications of disk features are pre- 
sented in Sect. ^ Section [S] focuses on the quantitative 
measurement of spiral features. The flattening distribu- 
tions of the disk features and galaxies are analyzed in 
Sect. The luminosity function and number fraction 
of dEs with and without disk features is the subject of 
Sect. 13 and the limitations in detecting disk features are 
considered in Sect.|Sl In Sect.|51we show how our objects 
are spatially distributed within the Virgo cluster, and a 
discussion and summary is given in Sect. ITni 

2. DATA AND SAMPLE SELECTION 

2.1. SDSS images 

The SDSS DR4 cov ers all galaxies listed in the Virgo 
Cluster Catalog fVCC. lBinggeh et alJ l985') with a decli- 
nation of (5 < 16?25, except for an approximately 2° x 2?5 
area at a « 186?2, S « -|-5?0 (see Fig.^l. It provides re- 
duced and calibrated images taken in the u, g, r, i, and z 
band with a pixel scale of 0'.'396, which corresponds to a 
physical size of 30 pc when adopting m — M — 31™0, i.e. 
d — 15.85 Mpc. The SDSS imaging camera takes data in 
drift-scanning mode nearly simultaneously in five photo- 
metric bands, u, g, r, i, and z, and thus combines very ho- 
mogeneous multicolour photometry with large area cov- 
erage, good resolution, and sufficient depth to enable a 
systematic analysis of early- type dwarfs. The images 
have an absolute astrometric accuracy of RMS < O'.'l 
per coordinate, and a relative accuracy between the r 
band and each o f the other bands of less than 0.1 pixels 
IjPier et alJl2003j) . They can thus easily be aligned using 
their astrometric calibration and need not be registered 
manually. The effective exposure time of 54 s leads for a 
bright dE (tob ~ 14) to a typical total signal-to-noise ra- 
tio (SNR) of about 1000 in the r-band within an aperture 
radius of approximately two half-light radii. For a faint 
dE (me ~ 18) this value is typically about 50. The RMS 
of the noise per pixel corresponds to a surface brightness 
of approximately 24.2 mag/arcsec^ in the u-band, 24.7 in 
g, 24.4 in r, 23.9 in i, and 22.4 in z. 

2.2. Image stacking 

In order to reach a higher SNR than that of the in- 
dividual images, we produced a combined image by co- 
adding the g, r, and i-band images. The u and z-band 
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images were not used, since their SNR is significantly 
lower and would thus lead to a decrease of the SNR of 
the combined image. When determining the sky level, 
proper object masks are required, so that pixels con- 
taining light from a star or a galaxy are excluded from 
the sky level calculation and only 'sky pixels' (i.e. pixels 
that contain nothing but sky background) remain un- 
masked. For this purpose, we applied the Source Extrac- 
tor Software (Bcrtin & Arnouts 199(|) to each object's 
image and each band to yield a 'segmentation image' 
which marks the pixels of all detected sources by assign- 
ing them non-zero values. To ensure proper masking of 
all objects, we expanded the source areas on the seg- 
mentation ima ge by smoot hing it with a Gaussian filter, 
using IRAF^ l)Todvlll993(l . The resulting image serves 
as object mask. The sky level was then determined with 
IRAF / imstat on the so-masked images along with the 
noise level, and was subtracted from the images. The 
g and i band images were shifted with IRAF / imshift to 
match the r band image; shifts were determined from the 
SDSS astrometry provided for each image (see above). 
We then applied wei ghts Wg^r.i to each image, following 
iKniazevet aLl(|20n4D : 
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with Sg^i-^i being the sky level and cTg.r.i the noise level. 
The weighted g, r, and i-band images were then summed 
to form the final combined image for each object. The 
resulting total SNR is about a factor of \/3 larger than 
in the r-band image. 

2.3. Sample selection 

From visual inspection of the combined images we 
chose a magnitude limit of jtib = 18'."0 for our study, 
with rriB provided by the VCC. This is the same magni- 
tude l imit up to which the VCC was found to be com- 
plete l|Binggeh et al.lll985j) . Adopting m - M = 3T.°0, 
it corresponds roughly to a limit in absolute magnitude 
of Mb < — 13™0. A more thorough examination of our 
limitations in detecting disk features is presented later in 
Section |S1 Initially, we selected all 552 cluster member 
and possible member galaxies with tub < 18'."0 that were 
classified as dwarf elliptical or dSO in the VCC, including 
those with uncertainties. We took into account the re- 
vised membership and classification from IBinggeli et al.l 
lfl993ljl . as well as updated cla ssifications for several ob- 
j ects given bv iBarazza et al.| (|2002, 2003), iGeha et aP 
l)2003t) . and lLotz et al.l 1)20041) . 25 galaxies are not cov- 
ered by the SDSS DR4. 25 objects with a classification 
'dE/dIrr' were excluded, and also all the remaining ob- 
jects were visually examined and excluded if they ap- 
peared to be possible dwarf irregulars due to asymmet- 
ric features in their image, which applied to 18 galax- 
ies. Thereby we avoid biasing our sample by the inclu- 
sion of potential non-early-type objects (which might be 
disk galaxies anyway). Three more objects (VCC0184, 
VCC0211, VCC1941) were classified as possible cluster 
members but appear to be probable background spirals 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



because of their small size and their spiral arm structure, 
and were therefore excluded as well. Five more objects 
(VCC0615, VCC0811, VCC1052, VCC1776, VCC1884) 
are of such low surface brightness that no examination 
for potential disk features is possible; these were also 
excluded. Our final sample comprises 476 early-type 
dwarfs, 414 of whic h are definite members of t he Virgo 
cluster according to lBinggeh et al.l l)1985lll993|) . 

3. IMAGE ANALYSIS TECHNIQUES 

For bringing to light weak features that are hidden 
by the dominating and mostly smooth and symmetric 
overall light distribution, two methods have proven suit- 
able. Unsharp masks are a common technique in detecing 
and enhancing weak substructure like e.g. nuclear bars 
or spirals (e.g. Liskcr et al. 2005; Erwin 2004). They are 
produced by first smoothing an image and then divid- 
ing the original by the smoothed one, which can easily 
be performed automatically on a large dataset. Another 
technique is to model the smooth axisymmetric light dis- 
tribution of a galaxy and subtract it from the original 
image ('e.g. IBarazza et al.ll200^ . with non- axisymmetric 
features like spiral arms remaining. Both methods 
have been used to identify spiral arms, bars, or edge- 
on disks in eight Virgo cluster early-type dwarfs so far 
(VCC0490, VCC0856, VCC094 0, VCCIOIO. VCC1036 
VCC1422, VCC1488 , yCC1695;LJ ericn ct al. 2000','200l 
I Barazza et al.l 120021: IGeha et all .2003: .Ferrarese et all 
l200fij) . These techniques are described below, along with 
the derivation of an elliptical aperture for each galaxy, 
which is required as input for both methods. 

From our ongoing analysis of dEs with blue central 
regions (paper II of this series) we know that a signif- 
icant fraction of dEs where no disk features were de- 
tected show obvious colour substructure. Since we ana- 
lyze the combined images from three bands in our search 
for disks, it could happen that colour substructure within 
the galaxy mimics the presence of a disk feature. To test 
this, we produced (uncalibrated) colour maps by divid- 
ing the aligned g and i-band images. Any detection of a 
disk feature with the methods outlined below can then 
be compared to the corresponding colour map and can 
thus be judged for reliability. To investigate whether or 
not there are any dEs in which colours do trace disk sub- 
structure requires a quantitative colour analysis that will 
be the subject of a future paper in this series. 

3.1. Elliptical apertures 

An elliptical aperture for each galaxy was determined 
by performing ellipse fits with IRAF / ellipse on the com- 
bined image, allowing center, position angle, and elliptic- 
ity to vary. One of the outer elliptical isophotes - usually 
between 1 and 2 half-light radii - was then chosen by eye 
to trace best the outer shape of each galaxy, as exempli- 
fied for VCC 1010 in the upper left panel of Fig. El This 
ellipse was adopted to define the ellipticity and position 
angle of the galaxy. 

3.2. Unsharp masks 

We produced a set of unsharp masks for each object by 
smoothing the combined image with a two-dimensional 
circular and elliptical Gaussian, one at a time, of var- 
ious kernel sizes a. A small value of a will enhance 
small structures and weaken large features at the same 
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time, while a large kernel size will enhance large struc- 
tures over small ones. For each set of unsharp masks 
we chose values of ct = 2,4,6,9,13,20, and 30 pixels. 
With d = 15.85 Mpc (m - M = 3r."0) and a subsequent 
pixel scale of 77 pc/arcsec (30 pc/pixel), these values cor- 
respond to 0.06,0.12,0.18,0.27,0.40,0.61, and 0.91 kpc, 
respectively. 

It is desirable to produce both masks created with 
a circular Gaussian (hereafter referred to as 'circular 
masks') and masks with an elliptical Gaussian ('ellipti- 
cal masks') corresponding to the galaxy's ellipticity and 
position angle. Circular masks of non-circular artifi- 
cial galaxies show a characteristic narrow shape along 
the major axis that could easily be confused with an 
edge-on disk and does not occur when applying ellipti- 
cal masks. We demonstrate this in Fig. |3| where a dE 
is represented by a two-dimensional exponential surface 
brightness profile with an elliptical shape created with 
IRAF / mkobjects (left panel). A circular unsharp mask 
with a Gaussian kernel of cr = 4pix, feigning an edge-on 
disk, is shown in the middle panel. In the right panel, an 
elliptical mask with position angle and ellipticity match- 
ing that of the galaxy has been applied: no substructure 
is seen. This is due to the fact that the scale radius of 
the light profile is smaller along the minor axis; therefore 
an isotropic Gaussian will blur the object much stronger 
along the minor than along the major axis. For detection 
of egde-on disk features or bars that are roughly paral- 
lel to the major axis, elliptical masks are thus clearly 
preferred. However, frequently the inner isophotes of an 
object are significantly rounder than the outer ones that 
define the Gaussian's ellipticity. In these cases, again 
an artificial narrow (bar-like) structure will appear along 
the minor axis, due to the very same effect as described 
above. Here, circular masks serve as a complementary 
check whether an apparent elongated feature along the 
minor axis is real or is only caused by varying ellipticity. 

3.3. Residual images from, ellipse fits 

A galaxy's surface brightness distribution can be mod- 
eled by performing ellipse fits (with IRAF / ellipse) and 
then feeding the output directly into the task bmodel. 
The resulting model image is then subtracted from the 
original object, yielding a residual image. Any infor- 
mation contained in the results of ellipse fitting directly 
enters the model. This can nicely be demonstrated on 
VCCIOIO, which hosts a bar. If we construct a model 
through ellipse fits with variable position angle and el- 
lipticity, the bar is not seen at all in the residual image 
(lower right panel of Fig. ^ since it has been fully re- 
produced by the model. If position angle and elliptic- 
ity are instead fixed at a value taken well outside the 
bar (namely the chosen elliptical aperture as described 
above), a strong residual double-cone is seen (lower left 
panel of Fig. O, which has already been explained by 
iBarazza ct al. ( 2002) as characteristic shape of a chang- 
ing position angle, and therefore of a bar. Similarly, spi- 
ral arms can be reproduced to a large extent by varying 
ellipses, and thus do not appear in the residual image 
unless position angle and ellipticity are kept fixed. 

From the above considerations it is obvious that any 
disk feature can best be detected with a model built 
through fixed ellipticity and position angle (later referred 
to as 'fixed model'). However, in principle any additional 



weak, asymmetric features would require variable ellipse 
parameters ('variable model'), so that the bar or spiral 
is properly reproduced in the model and fully subtracted 
from the image, and the additional substructure remains. 
Therefore, both types of residual images were visually ex- 
amined along with the unsharp masks for each object. 

3.4. Artificial galaxies 

In addition to the SDSS data we produced artifi- 
cial dE galaxies with IRAF / mkobjects, adopting a two- 
dimensional exponential surface brightness profile with 
an elliptical shape (left panel of Fig. I^jl. This 'primary' 
object was then superposed by another 'secondary' ex- 
ponential light distribution with the same or higher ellip- 
ticity, representing an (inclined) disk within a spheroid 
(Fig.0J. Various primary-to-secondary flux ratios, scale 
ratios, position angles and inclinations were reproduced, 
in order to provide a model counterpart for real galaxies 
that potentially are spheroids hosting a disk. The noise 
characteristics of the artificial images were chosen to be 
similar to a typical SDSS image, and galaxies covering a 
range of SNR values were created. 

4. RESULTS: EARLY-TYPE DWARFS WITH DISK 
FEATURES 

Close visual inspection of the combined image, the set 
of unsharp masks, and the two residual images was per- 
formed for each galax y, using the SAOImage DS9 tool 
Ij.Tove fc Mandel2003D . It turned out that unsharp masks 
are the primary means to search for substructure: espe- 
cially for small elongated features, they often provide 
a more reliable and clearer detection than the residual 
images do. In turn, only in very few cases did the resid- 
ual images show hints of substructure where the unsharp 
masks did not. However, in these cases the features 
were weak and their shape hard to define. Therefore we 
adopted a conservative approach and did not consider 
them as possible substructure. As IBarazza et al.l l)2003() 
pointed out, care must be taken with features seen solely 
on the residual images, since the models can be deceived 
by e.g. changing ellipticity and position angle, so that the 
resulting residual image would feign some substructure 
where none is present. Furthermore, the variable model 
turned out to be of little use, since it either reproduces 
substructure completely and yields a blank residual im- 
age (see Fig. |21), or leaves only weak features that are 
readily seen in the unsharp masks and the fixed model 
residual image. The situation described above that the 
variable model would bring to light secondary features 
by reproducing and subtracting the primary ones did 
not occur, i.e. no secondary substructure remained in 
the residual image other than weak and highly doubtful 
features. 

4.1. Disk detections 

We identified 14 out of 476 early-type dwarfs that un- 
ambiguously show disk features, as exemplified in the 
upper three panels of Fig. El Moreover, we find 'proba- 
ble disks' in 10 objects (third panel from bottom of Fig. 
13)), and 'possible disks' in 17 objects (lower two panels of 
Fig.|3J). This distinction between 'unambiguous', 'proba- 
ble', and 'possible' disks is based on the visual judgement 
of all three authors, and is intended to be an honest rep- 
resentation of the (un) ambiguity and the SNR of disk 
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features. In the case of a possible edge-on or inclined 
disk, we used comparisons with artificial two-component 
galaxies to check whether our interpretation is consis- 
tent with such a structure. This is exemplified in Fig. 
01 where the galaxy VCC0990 - classified as 'probable' 
dEdi - is compared to an artificial galaxy consisting of a 
'primary' and a 'secondary' component, the latter being 
fainter and having a larger ellipticity (i.e. representing a 
larger inclination angle). The simulated image is chosen 
to be similar in SNR and size, and indeed the shape of 
the galaxy images as well as their unsharp masks look 
similar. 

In two cases (VCC1684 and VCC1779), the colour 
maps (see Sect. EJ show a blue central region that is 
similar in appearance to the possible disk features. As 
a further test we produced unsharp masks for the two 
galaxies from the i-band images only. However, in both 
cases we can neither reject nor unambiguously confirm 
the presence of an inclined disk. We thus list both ob- 
jects as showing 'possible disk' features. 

In several cases we could not decide whether we see 
an edge-on disk or a bar; nevertheless, both were taken 
as disk feature, since the presence of a bar commonly 
requires a disk. Moreover, apart from the simple category 
'no substructure detected' (applying to 406 objects listed 
in Appendix Ie)| . we labelled 29 galaxies as objects where 
substructure of some kind is present, but not necessarily 
indicative of a disk ('other substructure'; objects listed in 
Appendix IE|| . 17 of these show irregular central features 
(also see Sect. 14 .211 . five have a boxy shape, in four objects 
a feature like a dust lane is seen, and for three objects 
the unsharp masks appear to show a luminosity excess 
in the inner part. 

Of the eight Virgo dEs for which disk features have 
been reported, five (VCC0490, VCC0856, VCCIOIO, 
VCC1036, VCC1695) are contained in our 14 unambigu- 
ous detections, and one (VCC1422) is a probable detec- 
tion. Both VCC0940 (reported by Barazza et all 12002(1 
and VCC1488 (reported by Gcha ct al. 2003) were not 
even identified as a dEdi candidate by us. The rea- 
son mi ght be twofold: first, those studies (as well as 
iDe Riicke et al . 2003) use a boxcar or median filter to 
create their unsharp masks. As we demonstrated above 
(see Fig. O and Sect. I3.2|l . applying such a filter to a 
perfectly smooth elliptical light distribution will yield an 
artificial elongated structure in the unsharp mask. This 
effect might well apply to VCC1488 with its axial ratio of 
0.55, but less likely to VCC0940 which has an axial ratio 
of 0.76. However, the disk features of both galaxies were 
also reported to be seen in the residual images result- 
ing from ellipse fits and subsequent modeling of the light 
distribution. Gi ven that bot h the data from Geh a et aD 
(|2003) and from lBarazza et al. ( , 2002') are of higher depth 
and resolution than our SDSS images, the non-detection 
of ours might simply reflect our limitations in detecting 
disks, and shows that more dEdis might exist than those 
identified by us (see also Sects. [3 and |H1). 

We list the dEdis and dEdi candidates in Table Q] We 
do, however, not attempt to reclassify objects, since clas- 
sification schemes in the VCC were fairly complex and 
based on the surface brightness distribution, whereas we 
aim solely at stating whether or not a dE's image shows 
features of a disk. In principle, it would be desirable 
to establish a 'pure' definition of the dSO class as those 



(and only those) dEs hosting (or being) a disk. Unfor- 
tunately, this is not possible: apart from the fact that 
many objects can only be termed candidates due to the 
limited SNR, those where no disk was found do not nec- 
essarily have to have no disk. It appears therefore most 
useful to not touch the original VCC classification, but 
instead to provide a list of (candidate) dEdis that can 
be correlated with all sorts of observables in future stud- 
ies of dEs. A thorough reclassification of all galaxies 
is deferred to a future study. We point out that our ob- 
jects are n ot related to the so-called dwarf spiral galaxies 
defined bv lSchombert et al.l l) 1995(1 : while those have a 
classical bulge, our objects do not. 

4.2. Correlation with the original dSO class 

iBinggeli fc CameronI l|1991l ^ described five cases in 
which a galaxy was classified dSO, with characteristics 
mostly indicative of a disk nature of the galaxy. Briefly, 
criteria for dSOs were a bulge-disk-like profile, high flat- 
tening, a lens-like appearance, a global asymmetry (like 
a bar or boxiness), and an irregularity in the central part. 

Our initial sample - prior to exclusion of possibly ir- 
regular objects - contained 47 out of 50 galaxies classi- 
fled as dSO or candidate dSO (e.g. 'dE or dSO') in the 
VCC. Two objects were then exluded due to a possi- 
ble irregular nature; thus 45 (candidate) dSOs are left in 
our working sample. 22 of these are indeed classified by 
us as dEdis or dEdi candidates, constituting 54% of our 
dEdi sample. 14 objec ts have 'other substructure' which 
reflects the criteria of IBinggeli fc Cam eron' ('1991'): 3 of 
them have a boxy shape, and 9 show irregular or clumpy 
central features likely caused by gas and dust. As an ex- 
ample for the latter, we show in Fig.Elthe image and un- 
sharp masks of VCC0781, which looks somewhat similar 
to the well-known dwarf elliptical NGC205 in the Local 
Group. Interestingly, all of these 9 objects with central 
gas/dust features have a blue central region with ongo- 
ing star formation or at least very young stars, similar to 
NGC205 and also to the galaxy presented bv iGu et alJ 
l|2006tl . This nicely confirms Binggeli's & Cameron's con- 
clusion, "the irregularity must stem from recent or ongo- 
ing star formation" (drawn without colour information 
or unsharp masks!). None of these galaxies shows (addi- 
tional) disk features; thus caution must be taken when 
treating them as dEdis only because of their dSO class: 
not all classified dSOs are dEdis. These objects might 
prove highly important for investigating possible forma- 
tion channels for dEs; therefore they will be the subject 
of paper II of this series (Lisker et al., in prep.). 

Finally, for 9 of the 45 (candidate) dSOs, neither a disk 
nor other substructure was found. However, three of 
these are classified 'dE or dSO', three are 'dSO?' (i.e. high 
uncertainty), and two are 'dSO:' (i.e. some uncertainty); 
hence we most probably did not miss any significant disk 
or irregular substructure. The one unambiguously clas- 
sified dSO (VCC1912) had been classified as such mainly 
due to high flattening. While our measured axial ratio 
of 0.33 is small, it is not small enough that we would 
classify it as dEdi based on flattening only. 

5. PROPERTIES OF SPIRAL FEATURES 

5.1. Relative strength 

For those three dEdis with the best-defined spiral arms, 
we now attempt to obtain an estimate of the relative 
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amount of light that constitutes the spiral arms, as com- 
pared to the smooth and axisymmetrically distributed 
light. We thus need to measure the flux of the resid- 
ual image (showing only the spiral arms) within a given 
aperture, and compare it to the total flux of the galaxy 
within the same aperture. We shall term this flux ra- 
tio the 'strength' of the spiral features. However, in the 
residual image the flux level in between the spiral arms is 
significantly negative: when fitting ellipses, the average 
flux value of each elliptical isophote is affected by the 
spiral arms and thus comes out slightly too high. Con- 
sequently, somewhat too much flux is assigned to the 
smoothly distributed light component, resulting in nega- 
tive flux values when subtracted from the original image. 
To avoid or at least minimize this effect, we obtain op- 
timized residual images through an iterative procedure 
outlined in detail in Appendix ^ yielding a lower and 
an upper limit for the strength of the residual features. 

The resulting residual images for our three dEdis are 
presented in Fig. [7| Note that it is not the case that our 
disk detections would have been more efficient if we had 
used such optimized residual images from the beginning: 
the contrast of residual features like spiral arms does not 
differ with respect to the initial residual images - only 
the average flux level is offset systematically. 

Apertures enclosing the spiral arms were now chosen 
manually, and the strength of the spirals was measured 
from the residual and the model flux within the same 
aperture. The nucleus as well as foreground stars or 
background objects were masked to avoid any bias. The 
results are listed in Table |21 VCC0490 has the strongest 
spiral features, which amount to 12-13% of the total light. 
The spiral of VCC0308 constitutes 8-11% of the light, 
and VCC0856 only reaches 6-8%. 

With these results at hand, we can now for the first 
time in the course of this paper address the question of 
whether dEdis are disk galaxies, i.e. are of flat oblate 
shape like VCC1304 (third row in Fig. [SJ, or whether 
they are spheroids hosting a disk component. The ratio 
of the light within the spiral features to the smoothly 
distributed light has been measured to be within 6-12% 
for our three galaxies. Therefore, when assuming that 
these objects are spheroidal galaxies hosting an embed- 
ded disk, the total light within the disk cannot be much 
larger than the light within the spiral features, since oth- 
erwise the disk would be the dominating component and 
the object would not be a spheroidal galaxy in the com- 
mon sense. Therefore, assuming the light within the spi- 
ral features to be of the same order as the total light of 
the disk component, the above ratio of 'spiral light' to the 
smoothly distributed light should be comparable to the 
ratio of the secondary to the primary component in our 
two-component model images. If, however, our galaxies 
would be genuine disk galaxies, the spiral features might 
well contain just a fraction of the total light of the disk. 
Consequently, if the disk is seen edge-on and compared 
to a suitable two-component model image, the ratio of its 
secondary to primary component should be significantly 
larger than the value measured for the (face-on) spiral 
features. Indeed, for those dEdis with apparent inclined 
disks that could not be confused with a bar, the sec- 
ondary component of the similar-looking model images 
is only 0.5-1 magnitudes fainter than the primary com- 
ponent, whereas the spirals measured above are 2.2-3.1 



magnitudes fainter than the smooth axisymmetric com- 
ponent. Although this is no final proof due to the small 
number of objects considered, it points towards dEdis 
being disk galaxies, instead of just having a disk com- 
ponent. Further arguments supporting this view will be 
presented in Sect. El 

5.2. A possible connection to faint SO/Sa galaxies 

Since the strengths of the three spirals measured above 
already differ within a factor of two, it might be inter- 
esting to see how the galaxies' images would appear if 
their spirals were stronger by a certain amount. For 
this purpose, we simply multiplied the residual images 
by a certain factor and added them to the model of the 
smooth component, thereby mimicking a stronger spi- 
ral. Strikingly, with only a 0.5-1 mag enhancement, the 
galaxy does not look like a dwarf elliptical or dwarf SO 
anymore, but instead like a spiral galaxy, although with- 
out a bulge. 

It might thus be no coincidence that more than a 
decade ago, one of us (B.B.) identified a handful of "faint, 
dwarfish looking SO/Sa" galaxies in the Virgo cluster 
(VCC0522, VCC1326, VCC1368, VCC1757, VCC1902) 
whose appearance is very similar to what has been just 
described (Fig.Cj). These objects differ from normal (i.e. 
giant) SO/Sa galaxies: their surface brightness profiles 
are similar to early-type dwarfs and remain flatter than 
the flattest possible King profile when going inwards, i.e. 
they apparently have no bulge (Binggeli, unpublished). 
Thus, they are hardly normal SO/Sa galaxies, which typi- 
cally have a high bulge-to-disk ratio. Instead they have a 
central luminosity excess just like the early-type dwarfs. 

One might thus term these objects 'dwarf-like SO/Sa' 
galaxies, to distinguish them from their giant counter- 
parts. A further investigation of their characteristics and 
a detailed comparison with early-type dwarfs will be the 
subject of a future paper in this series. For our present 
study, we selected those two with the best-defined spiral 
structure (VCC0522, classified Sa, and VCC1902, clas- 
sified SO/Sa), in order to measure the spiral strength 
like we did above and compare it to the dEdis. Their 
strengths turn out to be slightly larger than the aver- 
age value of the three dEdis and similar to the strongest 
dEdi spiral (VCC0490): 9%-13% for VCC1902 and 11- 
14% for VCC0522. Both objects are about half a magni- 
tude brighter than the brightest dE(di)s. It thus appears 
plausible that the dEdis and these objects belong to the 
same population of galaxies that extends to magnitudes 
brighter than those of dEs and differs from the 'classical' 
dwarf ellipticals. 

5.3. Pitch angle 

In order to confirm our above hypothesis, we measured 
the pitch angle of the spiral arms of both dEdis and the 
dwarf-like SO/Sa galaxies on the residual images. We 
used the method described by Ma (2001): a spiral arm 
is traced by manually selecting a series of image posi- 
tions that follow the arm. These are then fitted by a 
logarithmic spiral, taking into account the galaxy's in- 
clination and position angle which we adopt from our 
elliptical apertures (in the case of VCC1896 these val- 
ues were taken from the axial ratio measurement of the 
disk). We measured two arms of VCC0308, VCC0490, 
VCC0856, and VCC1896, one arm of each of the two 
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possible cluster members VCC0278 and VCC1671, and 
two arms of the two faint SO/Sa galaxies VCC0522 and 
VCC1902. The resulting values are shown in Fig. |Slas 
black symbols and compared to th e values for various 
Hubble types from iMa et all l|1999l^ (grey) . The dwarf- 
like Sa VCC0522 falls within the range of values of the 
dEdis while the dwarf-like SO/Sa VCC1902 Hes shghtly 
below. The dEdis best agree with Hubble type Sab/Sb 
while the dwarf-like SO/Sa galaxies - if taken together 
- fall in the range of type Sab. An independent check 
of our measurements is provided by iJerien et alJ l)2000D 
who find a pitch angle of 12? 1 for VCC0856. For the two 
arms, we derive the values 10?5 and 12? 1, respectively, 
thus being in good agreement with those VLT-data mea- 
surements. Our derived pitch angles are incompatible 
with spirals of very late type (>Sc), which are often con- 
sidered as potential progenitors for dEs; see Sect. 1101 for 
a discussion. 

6. FLATTENING DISTRIBUTION 

A flattening distribution for our galaxies can be ob- 
tained in two ways: for the disk features themselves by 
directly measuring or estimating their axial ratio, and 
for the galaxies as a whole, based on their ellipticities. 
The first distribution - which we shall term the flatten- 
ing distribution of the disks - serves as a basic test that 
the features we see are indeed disk features. This is of 
particular importance for the inclusion of 'probable' and 
'possible' disk features into our working sample of dEdis. 
In order to have a statistically significant sample, e.g. to 
derive the luminosity function (Sect. [Tj), we would like 
to include not only those dEs with unambiguous disk 
features, but also those with probable and possible disk 
features into our dEdi working sample. This requires the 
flattening distribution of disk features to be consistent 
with the assumption of an intrinsic flat oblate (and cir- 
cular) shape, which shall be examined in the following 
subsection. 

The flattening distribution of the galaxies - presented 
below in Sect. 16.21 - serves a different purpose: it will 
allow us to consider the question of the possible disk na- 
ture of the dEdis again. If they were spheroidal galaxies 
with a (weak) disk component, the distribution of axial 
ratios should be significantly different from that of disk 
galaxies. In turn, if their flattening distribution would 
be consistent with them having an intrinsic disk shape, 
they would very likely be genuine disk galaxies. 

6.1. Flattening distribution of the disks 

Although not possible with perfect accuracy, still an 
estimate of the inclinations of the disks (not the galax- 
ies) can be obtained from either the unsharp mask or 
the residual images. An ellipse was manually (by eye) 
fitted to the disk using that unsharp mask or residual im- 
age where the respective features stand out most promi- 
nently (exemplified in Fig.|51). The results are shown in 
the left panel of Fig. ^] as a running histogram (black 
lines), i.e. at each data point we consider the number 
of objects within the chosen bin- width of 0.1 (±0.05). 
We take into account all 36 dEdis and candidates that 
are certain cluster members. Galaxies where we cannot 
decide whether we see a bar or an edge-on disk were as- 
signed two values: a lower limit assuming an inclined 
disk (solid line), and an upper limit from the axial ratio 



of the galaxy as a whole, assuming the feature was a bar 
(dashed line). A theoretical distribution assuming a disk 
with an intrinsic axial ratio following a narrow Gaussian 
around a mean value fi — 0.25 with a = 0.01 and a 
randomly distributed inclination is shown as grey solid 
line for comparison (Mihalas fc Binnev 1981i'l . Within 
the expected uncertainties for our relatively crude mea- 
surements, the observed and theoretical curve are nicely 
consistent with each other. This strongly corroborates 
the hypothesis that the features we see are disks, and 
moreover, it supports the approach of including not only 
the unambiguous but also the candidate objects into our 
dEdi working sample for the purposes of our analysis. As 
a further test, we examined the flattening distribution for 
'possible disks' only - it turns out to be very similar to 
the distribution for all dEdis. It therefore seems plau- 
sible that most of our 'possible' disk detections actually 
are disks. Nevertheless we prefer to keep the term 'pos- 
sible' in order to reflect that uncertainties are present in 
our visual identification of disk features. 

6.2. Flattening distribution of the galaxies 

Based on the elliptical apertures described in Sect. 13. ll 
we put together the distribution of axial ratios of the 
(candidate) dEdis, shown in the right panel of Fig. 1101 as 
running histogram (black line). For comparison, we show 
the theoretical curve assuming an intrinsic axial ratio dis- 
tribution given by a narrow Gaussian with /i = 0.35 and 
(J = 0.02. Obviously, there is almost perfect agreement 
of observed and theoretical distribution, a compelling in- 
dication for an intrinsic disk nature of the dEdis! This 
view gains further support from the comparison with the 
distribution of dEs where no disk features were found 
(Fig. Ill|l : these objects are clearly consistent with a 
population of spheroids, and differ significantly from the 
dEdi distribution. It thus appears very likely that dEdis 
are genuine disk galaxies. A prototypical representation 
of how these disk galaxies appear when viewed egde-on 
might be given by VCC1304 (third row in Fig. O with 
its axial ratio of 0.3 2. ^_^ 

While Singgc li fc PopesciJ l|1995fl already found dSOs 
to be significantly flatter than dwarf ellipticals, the dif- 
ference is even more pronounced for our comparison of 
dEdis and dEs with no disk detection. This is explained 
by the fact that not all dSOs are dEdis and vice versa: at 
least some galaxies that were classified as dSO might be 
spheroids (see Sect. 14 .21) . 

The flattening distribution also allows us to test 
whether or not all bright dEs might actually be dEdis, 
but are not identified as such due to limitations of our de- 
tection method. When we modify Fig.llllsuch that only 
galaxies of the brightest one-(two-)magnitude-interval 
are considered (not shown), the distribution of dEs with 
no disk detection is inconsistent with all of them be- 
ing dEdis as well. Therefore, while we might miss some 
dEdis in our search for disk features as outlined in Sect. 
IHl we can exclude the possibility that all of the brightest 
dEs are disk galaxies - a significant number of objects 
need to be spheroids. 

7. DISK FRACTION VERSUS MAGNITUDE 

In the upper panel of Fig. El we show the distri- 
bution of dEs and (candidate) dEdis with respect to 
their B band magnitude provided by the VCC. For this 
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purpose we present our data as a running histogram 
with a bin- width of 1™0 (i.e. ±0™5). Only galaxies are 
considered that are certain cluster members according 
to Binggch ct al. (1985, 199^), resulting in 414 objects 
(light grey shaded histogram), containing 36 dEdis and 
candidates (dark grey shaded). The fraction of (can- 
didate) dEdis among all dEs is shown as black solid 
line: it reaches more than 50% for the brightest ob- 
jects, and then decreases to few percent at m^ > 16. 
Thi s 'disk fraction' might be of special interest, since 
e.g. lBinggeli fc CameronI 1)199 If) discuss a potential break 
in dwarf galaxy structure at Mb^t — ~16, which corre- 
sponds to niB = 15.7 given their m — M = 31.7. 

A plateau is seen in the running histogram (the lumi- 
nosity function) of our full dE sample, the position of 
which coincides very well with the location of the dEdis 
in the diagram. As a test, we subtracted the dEdi-counts 
from those of the full sample, but still a weak bump re- 
mains. However, we need to take into account the fact 
that we might have missed a significant number of disks 
in dEs due to the limitations of our data (which are as- 
sessed in Sect.|Hll. Therefore, we now multiplied the dEdi 
counts with 1.5 to account for the missed ones, and sub- 
tracted these counts from those of the full sample. In- 
deed, the plateau disappears (lower panel of Fig. IT^ . 

These results - independent of any considerations in 
previous sections - suggest very convincingly that dEdis 
are a different population than dEs with no disk, i.e. both 
have different origins not related to each other. Taken 
together with the indications for the disk nature of dEdis, 
evidence accumulates that dEdis are not just dwarf ellip- 
ticals with embedded disks, but instead constitute a pop- 
ulation of disk galaxies different and independent from 
classical dwarf ellipticals. 

In the following, we attempt to estimate the number 
of disks that are missed by our study, in order to assess 
whether the above assumption of a factor of 1.5 is real- 
istic. Moreover, we attempt to independently show that 
the decline of the disk fraction is real, and cannot be just 
an effect of limited data quality. 



8. LIMITATIONS IN DETECTING DISKS 

To obtain a realistic estimate for the limitations in 
detecting disks, we artificially dimmed our objects such 
that they correspond to dEs that are fainter by one (two) 
magnitude(s), also taking in to account the relation of dE 
magnitude and radius (.Binggeh fc Cameronlll99lh . This 
was done on the individual images and is described in 
more detail in Appendix IbI The resulting modified im- 
ages were then coadded like the original data, and un- 
sharp masks were created. The dimmed objects were 
then treated as if they were real galaxies that have to 
be searched for disk features, and the same categories 
('unambiguous', 'probable' etc.) were assigned. 

In Fig. 1131 we focus on the galaxies lying within the 
brightest one-magnitude-interval (solid line, filled cir- 
cles). When dimmed by 1 magnitude, they result in 
the histogram given by the dashed line with crosses, 
and when dimmed by two magnitudes, the resulting his- 
togram is shown by the dot-dashed line with filled tri- 
angles. This is illustrating the disk fraction we would 
expect to see at fainter magnitudes if the fraction of the 
brightest one-magnitude-interval of our sample would be 



constant with magnitude^. The obvious mismatch, along 
with the already strong decrease in disk fraction within 
the brightest one-magnitude-interval itself, suggests that 
the observed decline in disk fraction is real, and is not 
due to the limitations of the data. Even if we do not 
assume the true fraction to be constant, we find down to 
ruB ~ 16™0 the observed disk fraction to decline much 
stronger per one-magnitude-interval than what would be 
expected from artificial dimming (see Appendix [C)l . 

However, still a fair part of the decline is likely to be 
caused by the latter effect: the curve for objects dimmed 
by 1 magnitude lies at about a factor of 1.2 lower than 
the original one, and the two-magnitude curve is even a 
factor of two lower. This shows that our above estimate 
of the true number of dEdis being larger by 1.5 than what 
we observe is a useful estimate for the average fraction 
of missed objects. 

Still, the issue might be more subtle: if the relative 
strength of the disk features was decreasing with mag- 
nitude in addition to the SNR of the object as a whole, 
the estimate from artificially dimming the galaxies would 
be somewhat too high. While several of the disk fea- 
tures of the artificially dimmed galaxies would still be 
strong enough to be seen, some of the true observed ones 
would not. We examine this possibility in Appendix IdI 
and find that indeed somewhat more dEdis than esti- 
mated above might be missed at fainter magnitudes due 
to data limitations. However, if the true disk fraction 
were to decrease to zero this effect would be of minor 
relevance. Although we are not able to give an accurate 
estimate of the true number fraction of dEdis, we point 
out again that our analysis is consistent with the ap- 
proximation of multiplying the disk fraction with 1.5 in 
the lower panel of Fig. ^1 A significantly larger factor 
can be excluded following the argument given in Sect. 
16.21 the flattening distribution of the brightest one-(two- 
) magnitude-interval of our sample is inconsistent with 
all bright dEs being disk galaxies and instead requires a 
significant number of objects to be spheroids. 

9. SPATIAL DISTRIBUTION 

It is well known that the projected spatial distribu- 
tion of different morphological types of galaxies differs 
significantly ( the so-called morphology-density relation, 
"Dressler'1980) . Therefore it appears interesting to exam- 
ine the distribution of dEdis and dEs where no disk was 
found and compare it to other galaxy types. Those pro- 
jected spatial distributions are shown in Fig. ^^ along 
with the distributions for giant ellipticals (Es), for Es 
and giant SOs together, for spiral galaxies, and for irreg- 
ulars for comparison. Positions are taken from th e VCC 
by use of the VizieR database IjOchsenbein et alJ l20001 . 
Only certain cluster members are considered, and inter- 
mediate or uncertain classifications between the types are 
excluded^ Clearly, dEdis show the least clustering of all 
types, somewhat similar to the distribution of irregulars 
with 5 > 10°. 

For a more quantitative analysis, we show the cumu- 

■^ Here we neglect the fact that there is aheady a large decrease 
in disk fraction within the brightest one-magnitude-interval — how- 
ever, a certain interval width is necessary in order to still have a 
fair number of dEdis left among the two-magnitude-dinimed dEs. 

^ For example, a galaxy classified as 'E/SO' is excluded from the 
sample of Es, but included in the combined sample of Es and SOs. 
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lative distribution of each type of galaxy with respect to 
the distance from the cluster center. Since there is no 
unique definition for the latter, we decided to choose a 
point such that the radius of a circle enclosing all dEdis is 
minimized (Fig. ^ . For this purpose we use a 'corrected 
right ascension', which we define as 



(a - acenter) ' COs{S) + Ucenter 



(2) 



so that acorr is measured in true degrees. We choose 
'our' center to lie at 

Ctcenter = OLMS,7 — 0?15, Scenter = ^M87 — 0?85 (3) 

i.e. going from M87 slightly towards M86 and M49. In- 
terestingly, this circle at the same time encloses exactly 
all the giant ellipticals. For all other types, we only con- 
sider galaxies up to the maximum radius of the dEdis, 
in order to properly compare their clustering properties 
within that area. The fact that other galaxy types ex- 
tend slightly further outwards might have physical signif- 
icance, but could also be just due to the relatively small 
number of dEdis as compared to other types. Also, a part 
of it is due to the boundaries of the SDSS DR4 coverage, 
indicated with dashed grey lines in the upper panels of 

Fig.Hl 

We show the cumulative distributions in Fig.^] Along 
with the distribution for different morphological types, 
we show the expected distribution for an isothermal 
sphere {p{r) r^ r~^) in the upper panel, and for con- 
stant density (p(r) = const.) in the lower panel, where 
p{r) denotes the true volume density, not the projected 
surface density. This is done by populating a (three- 
dimensional) sphere at the distance of the Virgo cluster 
(taken to be d = 15.85 Mpc, i.e. m- M = 31'.°0) with 
the same number of objects as the number of dEdis, and 
then 'observing' the projected distribution of this sphere. 
Vertical intervals containing all but ±15.87% (= la) of 
Monte-Carlo-simulated values (darker grey) and all but 
±2.27% (= 2tT, lighter grey) are shown. Although a 
sphere is clearly not an ideal representation of the dy- 
namically young and unrelaxed Virgo cluster, this sim- 
ple model is intended to give at least a rough idea of the 
actual density distribution of the various galaxy classes. 

The well-known difference in the distribution of giant 
ellipticals and spirals or irregulars is clearly visible, and 
serves as guidance for the question of what constitutes 
a significant difference between two galaxy types in the 
diagram. The dEs where no disks were found roughly fol- 
low the distribution of giant Es and SOs, i.e. they are less 
centrally clustered than the Es alone but more strongly 
than spirals and irregulars. In contrast, dEdis lie clearly 
below this distribution, and for most of the sample they 
show even less clustering than spirals and irregulars, con- 
firming the impression given by Fig. 1141 While giant Es 
tend towards the isothermal sphere and spirals and ir- 
regulars more or less follow the distribution for constant 
volume density, dEdis even lie beyond the latter - a clear 
sign for them being not yet virialized, and thus being a 
population that has experienced fairly recent cluster in- 
fall. 

For the sake of completeness, we also show the result- 
ing distributions when M87 is chosen as cluster center in- 
stead (Fig. I16|l . The difference between dEdis and dEs is 
now slightly less pronounced, but also it now varies some- 
what less with radius than before. Here the dEdis follow 



closely the distribution of spirals, and fall within the 1- 
sigma area of the theoretical distribution for constant 
volume density - note, however, that a sphere around 
M87 is clearly no good representation of the Virgo clus- 
ter's shape. In contrast to the dEdis, dEs where no disk 
features were found approach the distribution of E+SOs, 
and at larger distances reach the distribution of the Es 
alone. 

10. DISCUSSION AND SUMMARY 

It is a long-standing question how early-type dwarf 
galaxies form, and whether there is more than one forma- 
tion channel producing them. Current theories include 
ram-pressure stripping, galaxy harassment, or in-situ for- 
mation. However, for a proper theoretical approach of dE 
formation, first the characteristics and possible subpop- 
ulations of the dE class need to be fully understood and 
unveiled from the observational side. While the defini- 
tion of the dSO class by BinEEcli ct al. (1985) implied a 
disk nature of these objects, the fairly diverse classifica- 
tion criteria had to remain suggestive but not compelling 
for dSOs being disk galaxies. The discovery of disk fea- 
tures in a handful of dEs had not yet been succeeded 
by a systematic, quantitative study, and could thus not 
provide significant input for models of dE formation. 

Moreover, kinematics - which might provide further 
insight into the presence of disks - are well studied only 
for a relatively small sample of dEs. With the SDSS data 
at hand, we performed for the first time a systematic 
search for disk features in an almost-complete sample 
of dEs down to me < 18™0, and found 41 out of 476 
objects to show (possible, probable) disk features. In 
the light of the diversity of the early-type dwarfs, one 
of our primary and most important results is that dEdis 
most likely constitute a different galaxy population than 
dEs where no disk features are found: the bump in the 
luminosity function of dEs (Fig. I12|l is highly unlikely to 
be an intrinsic characteristic of just a single population, 
and it is nicely explained by the superposition of dEdis 
and dEs with no disk features. Therefore, at least two 
different formation scenarios appear to be required: one 
for each dEs with and without disk features. 

When the first observations of spiral structure in dEs 
were made, galaxy harassment seemed to provide a sim- 
ple explanation for the apparently embedded disks in 
dwarf ellipticals: Mastrooictro ct al. ( 2005) showed that 
the progenitor galaxy's disk need not be completely de- 
stroyed during the process of transformation, but part 
of it is left over inside the newly formed dwarf. How- 
ever, we point out a main problem with this scenario: 
how could the observed well-defined, early-type spiral- 
arm structure of dEdis be reconciled with them having 
late-type progenitor spirals with their typically flocculent 
arm structure? Figure |S1 directly compares the pitch an- 
gle of our objects to that of Scd and Sd galaxies in the 
diagram - and shows an obvious mismatch. If one as- 
sumes a relatively weak spiral structure for the late-type 
progenitor that would quickly disappear after star forma- 
tion ceases, one might conclude that the above harass- 
ment scenario could still be valid, provided that the dEdi 
sp iral structure is pur ely of tidal origin, as e.g. suggested 
bv lJerien et al.l f2000) . The above question then changes 
into asking whether such well-defined spiral arms can at 
all be created through a process like harassment, and 
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what parameters determine their appearance. To con- 
firm that we are not looking at spiral structure traced 
by regions of star formation, we examined near-infrared 
H-band images for VCC0308 and VCC0856 which we ob- 
tained through the ESO/ST-ECF Science Archive facil- 
ity"* . These images show the very same spiral structure as 
the optical data, consistent with what would be expected 
for grand-design, early- type spiral arms. A detailed ex- 
amination of the colour properties of the spiral structure 
will be presented in paper III of this series (Lisker et al. , 
in prep.). 

Even without considering a specific formation theory, 
our data also allow us to address the question of whether 
dEdis are genuine disk galaxies or whether they are 
spheroids hosting a disk. Our distribution of axial ratios 
for the disks (i.e. where we measured the disk features, 
not the galaxies as a whole) agrees well with the expected 
distribution assuming an intrinsic axial ratio of 0.25 (left 
panel of Fig. I10|l . confirming our general approach to 
finding disk features in dEs. More importantly, also the 
distribution of axial ratios of the galaxies where disk fea- 
tures were found is nicely consistent with the assumption 
of their being disk galaxies with an intrinsic axial ratio 
of 0.35 (right panel of Fig. llU|) . This distribution signif- 
icantly differs from the distribution of dEs with no disk 
features, the latter being consistent with a distribution of 
genuine spheroids. We do not see how this could be rec- 
onciled with the assumption that dEdis themselves are 
spheroids - instead, we take these results as compelling 
indication for dEdis being disk galaxies, represented by 
the edge-on view of VCC1304 (third row of Fig.[SJ). 

Could this population of disk galaxies be simply an ex- 
tension of their giant counterparts? The deduced intrin- 
sic thickness of dEdis (0.35) agrees with t he correspond- 
ing value for giant Sa galaxies as given by Fouque et al. 
((1990) (0.37 for 'SO/Sa', 0.33 for 'Sa'; Schroder 1995), 
and the measured pitch angles best agree with Hubble 
type Sab/Sb (Fig.|S|). The dwarf- like SO/Sa galaxies pre- 
sented in Sect. l5.2l could in fact bridge the gap from dEdis 
to giant disk galaxies. Here, the presence or absence of a 
'classical' bulge can distinguish between what would be 
called a giant or a dwarf galaxy. However, early-type 
dwarfs are rare in the field environment, while early- 
type spiral galaxies are preferentially found in the field. 
This fundamental observation provides evidence against 
a close relation of dEdis and early-type spirals. 

The projected spatial distribution of dEdis within the 
Virgo cluster differs significantly from dEs with no disk 
features, and implies that the population of dEdis is not 
virialized yet. Thus, if dEdis would be the result of a 
morphological transformation, this should have occured 
recently. Any spiral structure of the late-type progeni- 
tor galaxies would have had to be destroyed during the 
process, since the spiral arm characteristics of the dEdis 
are incompatible with being remainders from Sc/Sd spi- 
ral galaxies. While a pure star formation origin of the 
spiral arms is unlikely (see above), they might originate 
from the recent galaxy-galaxy interaction that triggered 
the transformation process. Since such spiral structure 
would quickly disappear after the interaction ended, one 
would expect the dEdis to still show structural distor- 

* Observations made with ESO/NTT at the La Silla Observatory 
under programme ID 64.N-0288. 



tions, i.e. to be less homogeneous in appearance. More- 
over, a significant amount of tidal debris should still be 
present around them. At least the latter issue could be 
settled observationally with dedicated deep imaging of 
dEdis and their vicinity. 

Even before the discovery of the first spiral structure 
within a dE by iJerien et al.l l)2000() , it was obvious from 
the existence of a dSO class that treating all early-type 
dwarfs as one single population of galaxies always bore 
the risk of mixing objects that might have had differ- 
ent evolutionary histories. With our systematic search 
for disk features, we have now provided several strong 
indications that early-type dwarfs do indeed consist of 
two distinct populations of galaxies. Therefore, with our 
results at hand, we strongly recommend that those ob- 
jects identified by us as (candidate) dEdis be considered 
separately from the rest of dEs in any future study of 
early-type dwarfs, like e.g. a study of dE colours. Fur- 
thermore, one should keep in mind that a significant frac- 
tion of the brighter dEs where we did not find any disk 
features might still be dEdis - this possible incomplete- 
ness could fake systematic differences between brighter 
and fainter dEs. We also suggest to separately consider 
objects where we did not find disk features but that have 
been classified as dSO in the VCC, since our results con- 
firm that these also differ from 'ordinary' dwarf ellipti- 
cals. As a technical recommendation, we advise caution 
on the interpretation of substructure that is seen in un- 
sharp mask images created with isotropic smoothing of 
a non-circular object: as illustrated in Fig. El this can 
lead to artificial elongated features similar to an edge-on 
disk. 

Now that the separation between dEs and dEdis has 
been established, their properties can be analyzed. Given 
the disk nature of the dEdis, a correlation with kine- 
matical studies of early-type dwarfs is an obvious thing 
to do. Such a cor relation has first been investigated by 
iGeha et al.l l)2003f) who found that two out of three ro- 
tating dEs show disk features, yet two out of four non- 
rotating dEs have weak disk substructure as well. With- 
out going into the details of the kinematical analyses, we 
compiled results from several studies that state whether 
or not a dE show s signific ant rotation llvan Zee et alJ 
I2004H I(l^hr^ran f2003: Si mien fc Prugnie]ll2002l) . Note 
that these studies differ in their data properties, their 
maximum radius for sampling the rotation curve, and 
their criteria for significant rotation. 18 out of 29 galax- 
ies are found to be not rotationally supported, i.e. they 
show no or too slow rotation as compared to the ob- 
served velocity dispersion. 4 of these objects (22%) are 
(candidate) dEdis. However, it needs to be stressed that 
rotation curves are only sampled out to about the half- 
light radius, which might not be enough for definite state- 
ments about rotational support. 3 of those 4 dEdis do 
show significant rotation, but not enough to qualify for 
being rotationally supported. Of the 11 galaxies that 
were found to be rotationally supported, 6 (55%) are 
(candidate) dEdis. There is thus a tendency for dEdis to 
be rotationally supported systems, as one would expect 
for disk galaxies. The number statistics are consistent 
with our rough estimate of a third of the dEdis being 
missed in our study when assuming that most or all of 
them are rotationally supported. 

Given the different spatial distribution of dEdis and 
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dEs within the cluster, a further issue of interest would 
of course be their distribution of heliocentric velocities. 
These are available for 31 dEdis and 162 dEs where no 
disks were found. However, the two distributions do not 
differ significantly. Since the true threedimensional loca- 
tions of our galaxies within the cluster are not known, let 
alone the exact threedimensional structure of the cluster 
itself, unfortunately no useful conclusion can be drawn 
here. 

To demonstrate how our recommended separation of 
dEdis and the rest can be applied to other studies of 
early-type dwarfs, we show i n Fig. El the colours of dE 
nuclei derived bv iStrader et al.i (;200 5'l: five objects of 
this sample are identified by us as dEdis, and show redder 
nucleus colours than the bulk of dEs. To obtain a clearer 
relation, it would be desirable to further pin down the 
possible disk nature of the remaining dEs where we could 
not find disk features. This calls for a larger sample of 
kinematically studied dEs as well as for deeper images 
of higher resolution to detect further substructure, so 
that more quantitative input for theories of dE and dEdi 
formation can eventually be provided. 
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TABLE 1 
Early-type dwarfs with disk features. 



VCC 



me 



a.I2000 



Membership 



Note 



Certain disks 




1010 


13^72 


0523 


13.75 


2048 


13.85 


1036 


14.03 


0308 


14.30 


0490 


14.33 


0856 


14.42 


1695 


14.60 


1896 


14.78 


1671 


14.80 


0216 


14.90 


0278 


15.10 


1304 


15.50 


1204 


16.60 


Probable disks 




1422 


13.81 


1949 


14.19 


1947 


14.56 


1392 


14.62 


0407 


14.64 


0990 


14.81 


0218 


14.88 


2050 


15.20 


0336 


16.20 


1691 


17.30 


Possible disks 




1910 


14.17 


1183 


14.32 


0389 


14.40 


2019 


14.55 


0608 


14.70 


2042 


14.79 


1779 


14.83 


1684 


14.87 


1836 


14.92 


0397 


15.00 


1514 


15.10 


1444 


15.60 


0788 


15.80 


1921 


15.90 


2080 


16.20 


0854 


17.30 


1505 


18.00 



12^27™ 27?4 
12 22 04.1 
12 47 15.3 
12 27 41.2 
12 18 50.9 
12 21 38.8 
12 25 57.9 
12 36 54.9 
12 41 54.6 
12 36 32.2 
12 17 01.1 
12 18 14.4 
12 30 39.9 
12 29 38.0 



12 32 14.2 
12 42 57.8 
12 42 56.4 
12 31 55.9 
12 20 18.8 
12 27 16.9 
12 17 05.4 
12 47 20.6 
12 19 17.6 
12 36 51.1 



12 42 08.7 
12 29 22.5 
12 20 03.3 
12 45 20.4 
12 23 01.7 
12 46 38.2 
12 39 04.7 
12 36 39.4 
12 40 19.6 
12 20 12.2 
12 33 37.7 
12 32 35.9 
12 25 16.8 
12 42 26.5 
12 48 58.4 
12 25 55.7 
12 33 24.7 



+12°17'25" 
+12 47 15 
+10 12 13 
+12 18 57 
+07 51 43 
+15 44 42 
+10 03 14 
+12 31 12 
+09 35 05 
+06 10 11 
+09 24 27 
+06 36 14 
+15 07 47 
+07 06 24 



+10 15 06 
+12 17 14 
+03 40 36 
+12 10 28 
+09 32 44 
+16 01 28 
+12 17 22 
+12 09 59 
+05 52 33 
+12 57 31 



+11 45 15 
+11 26 02 
+14 57 42 
+13 41 34 
+15 54 20 
+09 18 27 
+14 43 52 
+11 06 07 
+14 42 55 
+06 37 24 
+07 52 17 
+09 53 11 
+11 36 19 
+11 44 25 
+10 35 12 
+12 46 11 
+15 24 28 



M 
M 
M 
M 
M 
M 
M 
M 
M 
P 
M 
P 
M 
M 



M 
M 
P 
M 
M 
M 
M 
M 
P 
M 



M 
M 
M 
M 
M 
M 
M 
M 
M 
P 
M 
M 
M 
M 
M 
M 
M 



3,4,(5) 

3,4,(5) 

1 

2 

5 

5 

5 

1,5 

3,5 

5 

5,(3) 

5 

2 

2 



1 

2,3,(4) 

3,4 

2 

2 

2 

2,(6) 

2 

1 

6,(5) 



1 
3 

4 

4,(5) 

2 

4,(5) 

2 

2,(7) 

5 

2,4,(3) 

2 

6 

2 

2 

2 

6 

6 



Note. — Objects are sorted by B-band magnitude ttib as given by Binggeli et al. (1985). Cluster membership is provided by 
IBinggeli et al J 119851119931) : M = certain cluster member, P = possible member. The last column contains information about the nature of 
the identified features; l = bar or edge-on disk, 2 = inclined disk, 3 = bar, 4 = disk, 5 = spiral arms, 6 = too flat for a spheroid, 7 = central 
gas or dust. The latter is an additional feature, but is not counted as disk. Numbers in brackets give uncertain features of which only a 
hint is present. 



TABLE 2 
Relative strength of spirals. 



VCC 


/res 
/mod 


Am 


/res 
/total 


(fc 


/ smoo 


Amsmoo 


\ /total / smoo 


0308 
0490 
0856 


0.107 
0.132 
0.075 


2.43 
2.20 
2.81 


0.097 
0.117 
0.070 


0.082 
0.122 
0.059 




2.71 
2.29 
3.07 


0.076 
0.108 
0.056 


0522 
1902 


0.159 
0.150 


2.00 
2.06 


0.137 
0.131 


0.127 
0.102 




2.24 

2.47 


0.113 
0.093 



Note. — Columns 2-4 give measured values for the optimized residual image without median smoothing, columns 5-7 give the same 
quantities for the version with smooting (see text for details). Columns 2 and 5 give the ratio of the flux of the residual image to the flux 
of the model image within the chosen aperture. Columns 3 and 6 give the same as a magnitude difference, and columns 4 and 7 give the 
fraction of residual to total light. 
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Fig. 1. — Distribution of dEdis within the cluster. Coordinates are given for J2000, and right ascension is corrected for the factor 
cos(<5), see text. Black circles are certain dEdis, black upward-pointing triangles are probable dEdis, and black downward-pointing triangles 
are possible dEdis. Grey crosses represent dEs where no disk was found. All other Virgo cluster galaxies with mB < 18™0 are shown as 
small black dots. Only certain cluster members are considered. The upper black cross gives the position of M87, the lower black cross 
marks our cluster center, chosen such that the radius of a circle enclosing all dEdis (dotted black line) is minimized (r = 5?0). Boundaries 
of the SDSS coverage are shown as grey dashed lines. 




Fig. 2. — Image analysis techniques. Upper left panel: Combined image of VCCIOIO, along with the elliptical isophote defining 
its shape. Upper right panel: Elliptical unsharp mask with kernel size cr = 20pix. Lower left panel: 'Fixed model' residual image, i.e. 
produced via ellipse fits with fixed ellipticity and position angle. Lower right panel: 'Variable model' residual image, i.e. produced via 
ellipse fits with variable ellipticity and position angle. Each panel has a horizontal scale of 300 pixels (119" or 9.13 kpc with d = 15.85 Mpc, 
i.e. m- M = 31™0). 




Fig. 3. — Circular and elliptical unsharp masks. Left panel: Simulated galaxy image created with IRAF / mkobjects, with exponential 
intensity profile, scale length along major axis 20 pixels, and axial ratio 0.5. Middle Panel: 'Circular' unsharp mask of the simulated galaxy, 
created with a circular Gaussian of kernel size cr = 4 pix. An elongated feature appears due to the application of a circular Gaussian to an 
elliptical object. Right Panel: 'Elliptical' unsharp mask with the same kernel size along the major axis, created with an elliptical Gaussian 
matching the position angle and axial ratio of the galaxy. Each panel has a horizontal scale of 138 pixels. 
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Fig. 4. — Simulated vs. observed dEdi. Upper panels: Combined image of VCC0990 along with its elliptical unsharp mask (cr = 9 pix). 
Lower panels: Simulated two-component galaxy image along with its elliptical unsharp mask (cr = 9 pix). The 'primary' component has an 
exponential intensity profile with scale length 30 pixels, axial ratio 1. The 'secondary' component has an exponential intensity profile with 
equal scale length, axial ratio 0.5, and a total magnitude 0™5 fainter than that of the primary component. The parameters are chosen to 
roughly match the appearance of VCC0990. Note that the simulation contains no nucleus, which is why the central region of the unsharp 
mask is brighter in the observed image than in the simulated one. Each panel has a horizontal scale of 248 pixels (98" or 7.55 kpc with 
d= 15.85 Mpc). 




Fig. 5. — Early-type dwrarfs with disk features: dEdis. Combined images and unsharp masks for three dEs with unambiguous disk 
features (top three rows), one probable dEdi (fourth row), and two possible dEdis (last two rows). The galaxies are, from top to bottom: 
VCC0308 (spiral arms; unsharp mask kernel size a = 20 pix), VCC1896 (bar and weak spiral arms; a = 13 pix), VCC1304 (edge-on disk; 
o- = 20 pix), VCC0990 (inclined disk, also see Fig.El a = 9 pix), VCC1183 (bar; a = 6 pix), VCC2019 (possibly inchned disk, maybe warped 
or distorted; cr = 13 pix). Each panel has a horizontal scale of 98" (7.55 kpc with d = 15.85 Mpc). 
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Fig. 6. — A dE with irregular central substructure. Combined image of VCC0781 (left panel), unsharp mask image with kernel 
size u = 4pix (middle panel), and unsharp mask with a = 9pix (right panel). Of those dEs where substructure other than disk features 
was found, this galaxy represents the subgroup of objects with central irregularities likely to be caused by gas and/or dust. Each panel has 
a horizontal scale of 46" (3.53 kpc with d = 15.85 Mpc). 




Fig. 7. — Residual images of spiral arms. Combined images as well as optimized residual images as described in Sect. 15. fl are shown 
for the three dEdis w ith t he best-defined spiral structure (VCC0308, VCC0490, and VCC0856 from top), as well as for the two dwarf-like 
SO/Sa galaxies (Sect.^J VCC0522 and VCC1902 (bottom). Each panel has a horizontal scale of 162" (12.48 kpc with d = 15.85 Mpc). 



16 



Lisker, Grebel, & Binggeli 



. ■ 1 1 1 ■ ■ 1 ■ , 1 1 ■ 1 1 ■ 1 ■ 1 i-i-m 1 1 1 

1 II 1 


1 
1 

1 


1 1 > . 

D 

K . 
Ill' 



25 - 



20 h 

D t5 
(L 10 - 



5 - 



SO/o So Sob Sb Sbc Sc Scd Sd dEdi dSO/o 
Type 

Fig. 8. — Pitc h angle vs. mo rphology. Grey symbols: Mean pitch angle and l-cr error bars of the spiral arms for various Hubble 
types as given by IMa et al.l 119991) . The value for type Sd was only derived from two objects. Black filled circles: Pitch angle for the 
certain cluster members VCC0308, VCC0490, VCC0856, VCC1896 (left to right) for two spiral arms each (connected symbol pairs). Black 
open circles: Pitch angle for the possible cluster members VCC0278 (left) and VCC1671 for one spiral arm each; in both cases the other 
arm could not be traced well enough. Black asterisks: Pitch angle for the dwarf-like SO/Sa galaxies (see text for details) VCC0522 and 
VCC1902 for two spiral arms each (connected symbol pairs; the values for the arms of VCC1902 are almost equal). 




Fig. 9. — Disk axial ratio measurement. Illustration of manual choice of a best-fitting elliptical aperture (right panels) for each disk 
feature. From top to bottom: VCC0490 (residual image), VCCIOIO (unsharp mask with kernel size a = 13pix), and VCC1304 (unsharp 
mask with a = 20pix). Each panel has a horizontal scale of 116" (8.95 kpc with d = 15.85 Mpc). 
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Fig. 10. — Axial ratio distribution for disks and galaxies. Running histogram with bin-width 0.1 for all 36 dEdis with certain 
cluster membership. Left panel: Distribution of axial ratio measurements of disk features as illustrated in Fig.|5] For the solid black line we 
assume that all elongated features where we could not decide between an inclined disk or a bar actually are an inclined disk. For the dashed 
black line, we assume that these features are bars, and thus adopt the axial ratio of the galaxy as an upper limit. The grey line shows the 
theoretical distribution for an intrinsic axial ratio represented by a narrow Gaussian of /.t = 0.25, a = 0.01, following iMihalas Hi BinnevI 
^JSl). It is normalized to the same area under the curve as the black solid line. Right panel: Distribution of axial ratios of the galaxies. 
The grey line represents an intrinsic axial ratio that follows a a narrow Gaussian of fi = 0.35, cr = 0.02, and is normalized like above. 
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Fig. 11. — Galaxy axial ratio distribution. Running histogram of the galaxy axial ratio distribution of dEdis (black) and dEs where 
no disk features were found (grey). Both histograms are normalized to an area of 1. 
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Fig. 12. — Luminosity function and disk fraction. Upper panel: Running histogram of the number of all dEs (light grey) and 
(candidate) dEdis (dark grey) with respect to B-band magnitude as given by the VCC. The bin-width is 1™0, therefore the counts are 
incomplete for m^ > IT^S (vertical dotted line). A bin is calculated at each position of a galaxy in the full sample. The upper x-axis 
gives absolute magnitude assuming m — M = 31™0. Only certain cluster members are considered. The ratio of both histograms is the disk 
fraction and is given as black symbols. Lower panel: Similar to the upper panel, but for all dEs minus 1.5 times the number of (candidate) 
dEdis (light grey), for 1.5 times the number of (candidate) dEdis (dark grey), and for the disk fraction resulting therefrom (black). 
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Fi g. 1 3. — Effect of SNR on the disk fraction. Upper panel: Running histogram of the disk fraction as given in the upper panel of 
Fig. 1121 fsohd and dashed Hne with circles). The brightest one-magnitude-interval is shown as solid line with filled circles, changing to a 
dotted line with open circles outside of the interval. When the dEdis in this interval are dimmed (see text) by 1 magnitude, the resulting 
disk fraction is given by the dashed line with crosses. A dimming by 2 magnitudes results in the dotted-dashed line with triangles. A 
histogram bin is calculated at each position of a galaxy in the full sample. 
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Fig. 14. — Distribution of morphological types within the cluster. For various types of galaxy (dEdis, dEs with no disk features, 
Es, Es-|-SOs, spirals, and irregulars) the projected spatial distribution is shown. Coordinates are given for J2000. Only certain cluster 
members are considered. The position of M87 is shown as grey cross. In the upper panels, boundaries of the SDSS coverage are shown as 
grey dashed lines. 
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Fig. 15. — Radial distribution of morphological types. Both panels show the cumulative distribution of the angular distances of 
galaxies from our chosen cluster center (see text). Only certain cluster members are considered, and all galaxy types are only considered 
up to the maximum distance of the dEdis. Various line types give the cumulative distributions for dEdis, dEs with no disk features, 
Es, Es+SOs, spirals, and irregulars, as labelled in the figure. In the upper panel, Monte-Carlo-simulations were performed to yield the 



expected distribution for an isothermal sphere potential (i.e. p 



'■') for a total number of 36 objects, i.e. the number of (candidate) 



dEdis. For the simulation, a distance to the Virgo cluster center of 15.85 Mpc was adopted (corresponding to m — M = 31), resulting in an 
angular scale of 0.28Mpc/°. The simulated objects populate a sphere with a physical radius of 1.4Mpc, i.e. corresponding to the angular 
value of 5?0 for the circle in Fig. The resulting distribution is shown as grey areas that enclose vertical intervals around the median, 
containing all but ±15.87% of simulated values (= Icr, darker grey) and all but ±2.27% (= 2(t, lighter grey). In the lower panel, analogous 
Monte-Carlo-simulations were done for a constant galaxy density. Note that l(2)-sigma areas are only valid for a comparison with the 
dEdis, not with other types, since the number of galaxies is different for the latter. 
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Fig. 16. — Radial distribution with respect to M87. Same as in Fig. 1151 but now adopting M87 as cluster center. 
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Nucleus colours vs. disk detection 
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Fig. 17. — Nucleus colours of dEdis. Histogram of g-z colours of dE nuclei (grey) as derived bv lStrader et all 120051) . Five of these 
objects are dEdis; their nucleus colours arc shown as black histogram. 
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Fig. A1. — Residual image optimization. Sketch of the iterative method for improving the spiral arm residual image. Each panel 
shows the flux distribution along an elliptical isophote, i.e. with respect to position angle. The isophote is represented by a constant flux 
value superposed by a crossing spiral arm modeled by a Gaussian. See text for the details of the method. Panel A shows the initial flux 
as solid line, with the average flux value given as grey line. Panel B shows the residual flux as dashed line, with the zero value given as 
dotted line. Panel C results from B when all negative values are set to zero. Panel D is obtained by subtracting panel C from A, with the 
new average flux value given as grey line. This value is subtracted from the original flux and results in the residual flux given in panel E. 
Panel F shows the flnal residual flux after 9 iterations. 



APPENDIX 
A. RESIDUAL IMAGE OPTIMIZATION 

In the residual images obtained in Sect. 13.31 the flux level in between the spiral arms is negative: when fitting 
ellipses, the average flux value of each elliptical isophote is affected by the spiral arms and thus comes out slightly 
too high (panel A of Fig. IA1|) . This results in negative flux values when the model is subtracted from the original 
image (panel B). We construct optimized residual images through the following iterative procedure. Where the initial 
residual image has negative flux values, its flux is set to zero, otherwise it is left unchanged (panel C). The resulting 
image is then subtracted from the original galaxy image (panel D), and a new residual image is obtained like before 
by fitting ellipses, constructing a new galaxy model, and subtracting it from the original image (panel E). This is 
repeated nine times iteratively, so that the final (tenth) residual image has reached (or come close to) a flux level of 
zero in between the spiral arms (panel F). A slight variation of this procedure is to smooth the residual image with a 
3x3 pixel median filter each time before the negative flux values are set to zero. It turns out that the final image of 
the latter version still has a slightly negative overall flux level, while the version without smoothing yields a slightly 
positive (i.e. too high) overall value in the residual image. We therefore use the strength measurement from the version 
with smoothing as lower limit, and the one without smoothing as upper limit. 

B. ARTIFICIAL DIMMING OF THE GALAXIES 

In order to artificially dim our objects by 1 (2) magnitudes, first the object size was decreased by a factor of 1.2 per 
mag nitude with IRAF / MAGN IFY, preserving the total flux. This follows the relation of magnitude and radius of the 
dEs ( Binggeh fc CameroiJll991|) : on average, the radius decreases with a roughly a factor of 1.2 per magnitude. Since 
this demagnification also affects the PSF, the image was then convolved with a (normalized) Moffat kernel of prope r 
size so as to approximately reproduce the original SDSS PSF (taken to be IFWHM = 4 pix: IStoughton et alJ l2002^. 
We then added noise to the image, with a a larger by 2.51 (6.56) compared to the original noise, thereby simulating 
the SNR of the 1 (2) magnitude fainter object. To increase the noise ct by 1 magnitude, one would actually need to add 
noise with a' = \j2.h\'^ — X'^-o — 2.30 o. However, since the original noise has already been weakened by demagnifying 
the image, we chose to use a' — 2.51 cr as a conservative approximation instead. 

C. EFFECT OF SNR ON THE DISK FRACTION 

In Fig. lCll we show a running histogram of the observed disk fraction (solid line with filled circles) and of the fraction 
obtained after dimming all objects by 1 magnitude (dashed line with crosses). The original disk fraction lies clearly 
below the shifted one until the region where both become very small and are affected by small number statistics. It is 
important to point out that Fig. ICll does not show how many dEdis would be detected assuming a constant true disk 
fraction. Instead, since all objects are dimmed by an equal amount (namely 1 magnitude), it shows the disk fraction 
that we would expect to find at a magnitude m when starting from the observed fraction at r?i — 1 and artificially 
dimming the objects there. Thus, any difference between the observed value at m — 1 and the 'expected' value at m is 
due to data limitations only. This is symbolized by the arrows in the figure. Consequently, if the observed decline from 
m— 1 to m is stronger than this 'expected' one, at least part of it has to be real and be not only due to data limitations. 
The ratio between the two curves thus tells us how much stronger the observed decline per one-magnitude-interval is 
compared to what artificial dimming of the galaxies would predict. We plot this ratio in Fig. ICll as grey solid line. 
Until ?7iB ~ 16™0 - where the number of dEdis becomes very small - for each one-magnitude-step the observed disk 
fraction declines a factor of 1.5-2.2 stronger than the 'expected' one from limitations of our data only. This clearly 
shows that the decline of the disk fraction is real. 
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Fig. CI. — Decline of the disk fraction. Running histogram of the disk fraction as given in the upper panel of Fig. 1121 fsolid line 
with circles). When all dEdis are dimmed by 1 magnitude (see text), the resulting disk fraction is given by the dashed line with crosses. 
A histogram bin is calculated at each position of a galaxy in the full sample. The grey line gives the ratio of both running histograms, 
and illustrates how much stronger the observed disk fraction declines per one-magnitude-interval than the 'expected' fraction does from 
artificial dimming only. 
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Fig. D1. — Effect of SNR on the flattening distribution. Upper panel: Distribution of axial ratio measurements of disk features as 
illustrated in Fig.|^ but this time only for the galaxies in the brightest one- magnitude-interval. The black lines show the running histogram 
of the original measurements, while the grey lines give the axial ratios measured for the disk features after artificially dimming the galaxies 
by 1 magnitude. A bin is calculated at each data point of each curve with a bin-width of 0.2, and the counts are normalized to an area of 
1 under the curve. Lower panel: Same as above, but here we compare the observed second-brightest one-magnitude-interval (black) to the 
brightest interval dimmed by 1 magnitude (grey, same as in the upper panel). 



24 



Lisker, Grebel, & Binggeli 



D. EFFECT OF SNR ON THE FLATTENING DISTRIBUTION OF THE DISKS 

Apparent axial ratios of disk features were not only measured on the original images, but also on those where 
the galaxies had been artificially dimmed, in order to reveal potential changes in the flattening distribution with 
magnitude. In the upper panel of Fig. IDll we compare the distribution of disk axial ratios for both the observed 
(black) and the artificially dimmed (grey) galaxies that lie within the brightest one-magnitude-interval of our sample. 
As in Fig. 1101 solid lines are derived from lower limits of the axial ratios, dashed lines from upper limits, depending on 
the interpretation of an elongated feature as a bar or as an edge-on disk. Both curves are normalized to an area of 1. 
They show a tendency for the dimmed objects towards lower axial ratios, indicating that the (artificial) dimming of 
objects might slightly prefer disks of certain inclinations over others. However, the distribution of axial ratios for the 
galaxies of the observed second-brightest one-magnitude-interval i s m uch more clearly skewed towards smaller axial 
ratios, i.e. larger inclinations (solid line in the lower panel of Fig. ID1() . While the upper panel suggests that part of 
this is due to the effect of the SNR on the detectability of features like spiral arms, it might also be that such features 
are intrinsically weaker - or not even present - in fainter objects. For example, the observed second-brightest one- 
magnitude-interval does not contain objects that look like the close-to-face-on spirals in all of VCC0308, VCC0490, 
and VCC0856, although it does contain galaxies with weaker spiral features that have a larger inclination. Note, 
however, that the black histogram in the upper panel consists of 16 objects, and both the grey histogram and the 
black histogram in the lower panel consist of only 13 objects. Therefore, the axial ratio distributions could at least 
to some extent be affected by small number statistics. We emphasize that the above effects on the axial ratio of the 
disk features need not go hand in hand with the axial ratios of the galaxies: as an example, the weak spiral arms in 
VCC1896 are not seen anymore when the galaxy is dimmed by one magnitude. One could then confuse the bar with 
being an inclined disk and thus measure a much smaller axial ratio of the disk feature, while the galaxy's axial ratio 
is the same in both cases. 



VCC numbers 
0011, 0029, 0033 



0200 
0361 
0503 
0652 
0761 
0824 
0928 
1039 
1111 
1198 
1296 
1386 
1453 
1539 
1647 
1711 
1815 
1895 
1971 
2063 



0208 
0372 
0504 
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0833 
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1040 
1115 
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0668, 
0765, 
0838, 
0931 
1044, 
1119 
1209 
1302 
1396 
1472, 
1553, 
1650, 
1719, 
1828, 
1901 
1982, 
2078 



E. OBJECTS WHERE NO DISK FEATURES WERE FOUND 

of objects where no substructure was found: 
0050, 0061 
0230, 0235 
0390, 0394 
0525, 0539 
0674, 0684 
0769, 0775 
0839, 0840 
0933, 0936 
1059, 1064 
1120, 1122 
1210, 1212 
1307, 1308 
1399, 1400 
1481, 1482 
1561, 1563 
1651, 1652 
1729, 1733 
1829, 1831 
1909, 1912 
1983, 1991 
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, 2081, 2083, 2088. 
VCC numbers of objects where substructure other than a disk 
0009, 0021, 0046, 0170, 0209, 0281, 0288, 0338, 0501, 0636, 0781, 
1457, 1501, 1512, 1567, 1617, 1668, 1715, 1743, 2045. 



was found (see Sect. l4.l|) : 

0870, 0929, 0951, 0962, 1078, 1288, 1334, 1370, 1395, 



